Anisotropic Coarsening: Grain Shapes and Nonuniversal Persistence 
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We solve a coarsening system with small but arbitrary anisotropic surface tension and interface 
mobility. The resulting size-dependent growth shapes are significantly different from equilibrium 
microcrystallites, and have a distribution of grain sizes different from isotropic theories. As an 
apphcation of our results, we show that the persistence decay exponent depends on anisotropy and 
hence is nonuniversal. 
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The geometrical Wulff construction [Q gives an explicit 
relation between the anisotropic surface tension and the 
resulting equilibrium crystal shape. This marks an early 
and dramatic success in quantitatively connecting mor- 
phology to the interfacial properties of a material. How- 
ever, distinct WulfF microcrystallites must be in 'splen- 
did isolation' — with negligible exchange between them 
in comparison to the internal dynamics required to equi- 
librate In contrast, dilute phase separating alloys 
and coarsening polycrystallitcs exhibit growing micro- 
crystalline droplets or grains with non-negligible inter- 
actions. While it has been shown for these and other 
coarsening systems that anisotropy influences the mor- 
phology 1^,0, such eflFects have not been quantitatively 
understood for even the simplest models of curvature- 
driven growth. 

The understanding of interacting isotropic phases (see, 
e.g., [||) was significantly advanced by the models of Lif- 
shitz and Slyozov 1^ and Wagner for diffusive and 
curvature-driven coarsening, respectively. These mean- 
field theories correctly capture a remarkable amount 
of coarsening phenomenology, and are exact in the di- 
lute limit. With this inspiration, we generalize Wag- 
ner's model — an interacting ensemble of coarsening 
droplets, evolving to continually lower their surface en- 
ergy without changing their total volume — to include 
arbitrary anisotropy in the surface tension and the in- 
terface mobility. We solve the model perturbatively in 
anisotropy strength, and relate the interfacial properties 
to the resulting non-trivial grain shapes. These "growth 
shapes" are contrasted with those of equilibrium (Wulff- 
constructed) grains to highlight the connection between 
dynamics and microcrystallite morphology. We then 
compare our results on the ensemble of grains to Wag- 
ner's isotropic solution to demonstrate anisotropy effects 
on coarsening correlations, including the effect on persis- 
tence exponents. 

Our model is applicable to single-phase polycrystallite 
coarsening, where distinct grains are distinguished only 



by their crystallographic orientation (see ||8|-[l0t|). Most 
theoretical studies of polycrystallites focus on their cellu- 
lar structure, specifically on the static and dynamical de- 
scription of the vertices where three or more grain bound- 
aries meet. However, vertex-based models have signifi- 
cant shortcomings when anisotropy is included, since it 
modifies both the distribution of the number of vertices 
per grain and the otherwise fixed angles formed where 
three grains adjoin Furthermore, von Neumann's law, 
a direct relationship in two dimensions (2D) between the 
number of vertices per grain and its area growth rate , 
no longer applies. With anisotropy, the evolution of a 
grain's area requires the complete specification of grain 
shape — including the orientations and, in general, the 
non-uniform curvatures of the interfaces. 

We present a complementary vertex-free approach to 
examine grain shape via an anisotropic dynamical mean- 
field theory. The neighboring grains outside the grain 
of interest are treated as providing an isotropic mean- 
field. We retain the crystallinity of the grain through an 
anisotropic surface tension and interface mobility, which 
results in the anisotropic Wagner theory. (A similar con- 
nection can be made between isotropic Wagner theory 
and soap froths |11|.) Ultimately, a synthesis of the 
present work with vertex-based models is desirable [ p^ . 

We find a dynamical scaling solution typical of coars- 
ening systems including clean polycrystallites. The 
characteristic length scale grows as a power law, L ^ i^/^, 
as expected for curvature driven growth jl^ . In the scal- 
ing regime the initial conditions are "forgotten," and the 
morphology, when scaled by the growing length, L(t), 
is invariant. Grain shapes of particular scaled size are 
also time- independent. These growth shapes are gener- 
ally quite different from equilibrium Wulff shapes — even 
when the mobility is isotropic! The isotropic grain size 
distribution is also modified by anisotropy, as discussed 
later. 

With our results, we can answer the question of uni- 
versality in persistence decay exponents. The persistence 
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is the fraction of the system that has not been crossed 
by a domain wall up to time t |T^-p^ . The decay of per- 
sistence to zero, P ^ , even from a starting time deep 
within the scaling regime, implies that every point in the 
system will eventually "realize" that equilibrium has not 
yet been reached. Persistence decay is a local signature of 
the non-equilibrium dynamics of the system. The degree 
of universality of this dynamical exponent has remained 
an open issue since no precise results have been obtained 
for models with non-trivial temperature dependence [ p^ . 
(Simulations have not yet found any temperature depen- 
dence within their accuracy (l^.) Since anisotropy varies 
with temperature, our model provides such a non-trivial 
temperature dependence in a coarsening system that we 
can then analytically relate to the resulting structure and 
to the persistence exponent, We find that 9 de- 

pends on both the anisotropics of the surface tension and 
of the interface mobility, so the persistence exponent is 
nonuniversal in anisotropic systems fl^ ]. 

We restrict ourselves to 2D, where the surface tension 
cr{ip) and the interface mobility M{ip) may be defined 
in terms of the angle ip between the interface normal 
and an arbitrary crystallographic axis. The anisotropic 
Allen-Cahn equation 20| is then derived from the lin- 
ear response of the interface to the local drive given by 
the Gibbs-Thompson condition, [a{ip) + a" {ip)]K, were a 
is the surface tension, and k the local interface curvature. 
The stiffness, a + a" , reflects the local change of extent 
and orientation of the interface due to a deformation. By 
allowing the interface mobility to depend on orientation 
pl| , and by including an applied field A coupled to one 
of the phases, we obtain the normal interface velocity 
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servation of the total area of all of the grains uniquely de- 
termines X{t), resulting in precisely the anisotropic Wag- 
ner theory. 

To proceed, we Fourier expand the anisotropic surface 
tension and mobility. 



a{tp) = (To 1 + <5 ^{o-fc cos{kifj) + dk sin(fc?/')} 



(2) 



k=l 



M{^) = Mo 



1 + (5 ^^{rrifc cos(fc?/;) -I- mfc sin(fc'0)} , (3) 



fc=i 



where S is introduced to organize a perturbative cal- 
culation. Wc parameterize each grain by a polar ra- 
dius R{(t>), as depicted in Fig. |l|, from which the in- 
terface orientation follows: 'ipi't') = 4' ~ arctan(i?'/i?) 
where R' = dR/d<j). Considering only smooth grain 
profiles, we relate normal and radial growth velocities, 
Vr{4>) = Vny/T+(~Rf/R)^, and calculate the curvature 
= [i?2 + 2i?'2 - RR"]/{R^ + i?'2)3/2. We then ex- 
pand i?. 



i?(0) = i?o 1 + ^{pk cos{k(j)) + pk sm{k(j})} 
with coefficients 



k=0 



Pkix) = akix)S + bkix)S^ 



(4) 



(5) 



and similarly for pk p^ . Grain sizes are labeled with 
a reduced length x = Ro/L, where L = {Moaot/2y/^ . 
For S — we recover Wagner's isotropic theory, with the 
familiar distribution of grain sizes (see [^jl^jlj]): 



fix) = ei^2xexph4/(2 - x)]/(2 - x)' 



(6) 



FIG. 




1. An anisotropic drop illustrating 



the polar co- 
ordinate, and V'(0)i the angle of the interface normal at the 
point (^R{(f>), 0). 

We now consider an ensemble of polycrystallite grains. 
Our mean-field approximation entails keeping only the 
crystalline anisotropy of each grain (ignoring its neigh- 
bors), neglecting vertices, and determining a self- 
consistent mean-field A to represent the effects of neigh- 
boring grains that may be growing or shrinking. The con- 



(The e prefactor is the area fraction of a randomly se- 
lected subset of grains — used later to calculate per- 
sistence.) For convenience, we define Rq by the re- 
quirement that X maintains this isotropic grain-size dis- 
tribution up to an anisotropy-dependent normalization. 



F2 = ^2°"* -I- S'^F2^' + — This requirement leads to non- 
zero po terms in the expansion (^ but preserves the range 
of scaled sizes, x £ [0,2]. (Note that po = 0.) Physical 
length scales, such as the grain perimeter, can be consis- 
tently derived from our results, as discussed below. 

The resulting interface equations for the ensemble of 
grains may be solved order by order in S [ p3| . The ze- 
roth order results reproduce the isotropic theory; the first 
order equations are new, and serve to determine a size- 
dependent grain shape through ak{x): 



a 17(2) 



x{2 - xfa'kix) - 4(A:2 + x- 2)ak = 
- k^)ak + '^{l - x)mk, 



(7) 



for fc > 1, with an identical equation for aj, in terms of 
(Tfc and rfik- For k > 1 the solution is 
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ak{x) =mk + K -mfc)(l - l/k'^)[l + flik^v)] (8) 

where n{k, v) = 2r(2 - k'^ ,v)v^^ , and v = k^x/{2 - 
x). We also have ao{x) = ai{x) — 0, the latter by our 
choice of coordinate origin . Clearly the grain shapes 
depend on grain size, through fJ. Even when the sur- 
face tension is isotropic {ak — for all fc > 0), we can 
obtain anisotropic grain shapes through the interface mo- 
bility. This is illustrated in Fig. ^ for a particular choice 
of M(V'). 
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FIG. 2. First order grain shapes for various sizes (not to 
scale) with an isotropic surface tension (tq but with a par- 
ticular anisotropic mobility 5m2 = 0.4 and 5rrn — 0.9 (all 
other rrifc — rhk = 0). The scaled grain sizes are, from the 
innermost, x = 0.01,0.5, 1, 1.5, and 1.99. For no value of x 
is there a circular grain, the equilibrium Wulff shape. Angles 
for which M{il)) is larger correspond roughly to larger radius 
in growing (larger) grains, and smaller radius in shrinking 
(smaller) grains. 

At all orders of S the equations for the grain shape are 
similar to (^, although the right-hand side will include 
products of lower-order solutions. While these equations 
are progressively more difficult to solve, we can itera- 
tively demonstrate that the solutions are finite at every 
order of 5 

In the special case where ruk — Uk for all k grains of all 
sizes have the equilibrium Wulff shape. This result holds 
to all orders in (5, and is due to a remarkable symmetry 
held by the interface equation (^ . The equilibrium grain 
shape is given by 



i?cq(</') 



cos 



(9) 



For this Wulff shape, a variational calculation shows that 
+ ct" {tp)]K is independent of angle [p3|| , from which 
we obtain Vr oc M(-0)-\/l + {R' / R)^ for all angles. If and 
only if the dynamical mobility anisotropy equals that of 
the static surface tension — that is, M{'ip) oc ct(V') — 
then we recover Vr oc i?oq ((/>), the condition for Wulff 
grains to keep their shape while evolving. This symme- 
try, evident in (^, leads to size- independent drop shapes 
and also shows up in the drop size distribution and per- 
sistence results, as discussed below. 

However, the dynamic mobility and the static surface 
tension will not be proportional except by special con- 



struction. Regardless, in physical systems M and a have 
different temperature dependences so that equality could 
not be maintained as temperature varies. In the general 
case, we will have size-dependent drop shapes given, to 
first order, by (^). In comparison, the Wulff construc- 
tion gives = ak as the leading contribution to the 
equilibrium grain shape. Even with an isotropic mobil- 
ity, mfc — 0, growth shapes differ from equilibrium and 
depend on grain size. 

The isotropic grain size distribution (^ applies only 
to our index Rq ||2^. Physically relevant lengths, such as 
extracted from the grain perimeter or the area, will gen- 
erally have different distributions. For example, the area 
A = i d4>R{(t)Y can be used to define Ra — \/ A/tt 
where 



Ra — 



(10) 



A scaled size z — Ra/ L may then be introduced, which 
wiU be related io x hy z ^ x + 6'^h{x) + 0{S^). The 
"area radius" distribution g(z) is then determined by 
g{z)dz = f{x)dx so that 

g{z) = f{z) - 5'[f{z)h'{z) + f'{z)h{z)] + 0{5') (11) 

where f{x) is the isotropic distribution p^ . The grain 
perimeter distribution follows similarly, though with a 
different function h{x). [In the special symmetric case, 
where ruk — (Jk for all fc > 1, all physical lengths have 
the same distribution. Since the grain shapes are size 
independent, h{x) — Jiq, and g{z) and f{x) differ only by 
an overall normalization.] 

We may also calculate the slow decay of persistence 
due to the evolution of a small area fraction e of randomly 
chosen grains, following [Him. The persistence P> of 
the region outside the chosen grains decays due to grow- 
ing grains via 9tP> = — i'>P>. The rate of encroachment 
of growing grains, w>, can be calculated from the grain 
shapes and (0). The power-law decay of persistence fol- 
lows directly from the result w> oc 1/t, with persistence 
exponent 9 — tvy. Anisotropy appears at 0{S^). The 
calculation is lengthy and details are reported elsewhere 
[p3|; however, the result simplifies to 



= "70 



(12) 



where 6*0 — 0.48797e is the 2D persistence exponent for 
the isotropic case [0 . We find that 9 equals the isotropic 
value 9o only when M{iIj) oc (t(-!/') (this holds to all orders 
due the symmetry mentioned earlier) and differs from 9q 
for any other anisotropic conditions. The order e coefR- 

(2) 

cients 0^, may be determined by numerical integration, 
and are well approximated by a large k expansion [p5[. 
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The persistence exponent depends continuously on both 
the mobihty and surface tension, and consequently on 
the temperature. The 2D Ising model provides an ex- 
plicit example: the anisotropic surface tension is known 
analytically for < T < Tc and the anisotropic mo- 
bility is known close to T = for Glauber dynamics p7| . 
At T = 0, and using only the leading contribution (12), 
we find 9 ~ 1.0344 6*0. The effect on the exponent is 
small but non-zero, and may be detectable numerically 
with more accurate studies. 

In conclusion, we have constructed a mean-field model 
for 2D polycrystallite coarsening with anisotropic surface 
tension and mobility. We find an exact scaling solution 
with size-dependent grain shapes that are generally unre- 
lated to the equilibrium Wulff shape. We use our solution 
to calculate the exponent describing persistence decay, 
and find that it is continuously dependent on anisotropy 
and hence nonuniversal with respect to temperature . 
We expect similar results to hold in three-dimensional 
systems. 

We hope that this study stimulates further research of 
the connections between nonequilibrium structure and 
anisotropic mobility and surface tension. Our next 
step will be to develop the anisotropic generalization of 
Lifshitz-Slyozov diffusive coarsening in bulk systems 
We also feel the influence of anisotropy on nonequilib- 
rium exponents needs further study. For example, we 
suspect that persistence exponents and their various gen- 
eralizations (see, e.g., ||l^) will prove to be nonuniversal 
whenever correlation functions are anisotropy dependent. 
Finally, we stress that persistence decay still provides an 
important description of nonequilibrium dynamics, and 
remains universal in intrinsically isotropic systems such 
as binary fluids, polymer blends, and soap froths. Per- 
sistence is particularly useful in discriminating between 
different dynamical models and in probing the dynamics 
of soap froths (l|j2|. 
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